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Motivated by recent experiments on large quantum dots, we consider the 
energy spectrum in a system consisting of N particles distributed among 
K < N independent sub-systems, such that the energy of each sub-system 
is a quadratic function of the number of particles residing on it. On a large 
scale, the ground state energy E(N) of such a system grows quadratically 
with N, but in general there is no simple relation such as E(N) = aN + bN 2 . 
The deviation of E(N) from exact quadratic behavior implies that its second 
difference (the inverse compressibility) xn = E(N + 1) — 2E(N) + E(N — 1) 
is a fluctuating quantity. Regarding the numbers xn as values assumed by a 
certain random variable Xi we obtain a closed- form expression for its distri- 
bution F{x)- Its main feature is that the corresponding density P(x) = dF }^ 
has a maximum at the point x = 0- As K — ► oo the density is Poissonian, 
namely, P{x) e~ x . 
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I. MOTIVATION 



Statistics of spectra is an efficient tool for elucidating properties of various physical 
systems. So far, most of the effort is focused on the study of energy levels of a system with 
a fixed number of particles. In this context, one of the central earlier results is that the 
spectral statistics of many-body systems such as complex nuclei agree with the predictions 
of random matrix theory On the other extreme, it was found that level statistics of a 
single particle in chaotic or disordered system also obeys a Wigner-Dyson statistics 0,§J. 

Recently, experiments are designed to get information on the statistics of the addition 
spectra of electrons in quantum dots ||. The pertinent energy levels E(N) are the ground 
state energies of a system consisting of N electrons residing on a quantum dot, which is 
coupled capacitively to its environment. 

Let us single out two properties of the addition spectra of quantum dots. The first one 
is that, on a large scale, the energy E(N) grows quadratically with N, while the second one 
is a consequence of charge quantization, namely, there is, in general, no simple relation such 
as E(N) = aN + bN 2 . In this context, an appropriate quantity whose statistics is of interest 
is then the inverse compressibility, 

Xn = E(N + 1)-2E(N) + E(N-1). (1) 

It is the deviation of E(N) from exact quadratic behavior which makes its second difference 
Xn non-constant. Indeed, in a recent experiment on large quantum dot it was found that 
the inverse compressibility vanishes for numerous values of electron number N. 

In the present work we study the statistics of the addition spectrum of a simple physical 
system with the two basic properties mentioned above. One example of such a system is 
motivated by considering the electrostatic energy of large quantum dots (although it should 
be mentioned that the model is too simple to describe the actual physics). To be specific, 
we have in mind a system of K metallic grains such that the number of electrons on the i th 
grain is rii {i — 0, 1, 2, . . . , K — 1) and their sum equals N. The electrostatic energy of the 
pertinent system is a bilinear form in the numbers rii with a K x K matrix w = \C~ l . Here 
C is a positive-definite symmetric matrix of capacitance and inductance coefficients. If the 
metallic grains are very far apart, the matrix C is nearly diagonal. Thus, we concentrate on 
the special case C = diag[Cj], for which the energy of the system is given by 

A'-l ^ 

E(N) = min £ — n?, (subject to Elt n i = N). (2) 

i=0 ALy i 

The minimum in (0) is taken over all possible partitions rii of N. 

Another example is the energy of a system composed of K different harmonic oscillators, 
among which one distributes N spinless fermions. If there are fermions on oscillator 
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i (whose frequency is u>i), then the energy of this oscillator (up to a constant) is Ei = 
TiujiTiiirii + 1), and hence the ground state energy of the system is 

K-l 

E(N) = min ^ E h (subject to Efjo 1 n i = N )- ( 3 ) 

i=0 

We will concentrate on the first example, which is borrowed from the electrostatics of quan- 
tum dots @, and refer to the constants Cj as capacitors. Some results pertaining to the 
second example (the system of oscillators @) are also presented. 

Regarding the numbers xn of (jl]) as values assumed by a certain random variable, the 
distribution of this random variable is the main focus of the present work, which culminates 
in Theorem 1, where we find a closed-form expression for the distribution. 

The problem of elucidating the (addition) spectral statistics of a a many-body system, 
consisting of several independent sub-systems (whose dependence of E on n« is known), 
looks deceptively simple. As will be evident shortly, this is not the case, and finding the 
distribution in question is quite a non-trivial task. Note that, even for a single particle 
system composed of several independent sub-systems (e.g., a system of a particle in several 
boxes), the derivation of level statistics requires a large degree of mathematical effort [|7|]. 
The rest of the paper is therefore devoted to a rigorous derivation of our main results. 



II. FORMALISM 

Definition 1. Let (^ n )^i be a sequence of real numbers and F a distribution function. 
The sequence (9 n ) is asymptotically F -distributed if 

\{l<n<M :d n <x}\ 

*• r (X) 

for every continuity point x of F (where \S\ denotes the cardinality of a finite set S). 

An equivalent condition is the following. Denote by St the point mass at t, and let \x be 
the probability measure corresponding to the distribution F (namely, n(A) = J l^dF(x) for 
any Borel set A). Then (8 n ) is asymptotically F-distributed if 

-(Se^Se^.-. + SeJ^fM 

(the convergence being in the weak*-topology). 

The notion of asymptotic distribution has a stronger version whereby, instead of requiring 
only that initial pieces of the sequence behave in a certain way, we require this to happen 
for any large finite portion of the sequence. This leads to 
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D(A) = lim 



Definition 2. In the setup of Definition 1, (9 n ) is asymptotically well F- distributed if 

\{L<n<M :9 n <x}\ . 
— * t ( X J 

M - L M-L^oo V ' 

for every continuity point x of F. 

Recall that the density of a set A C N is given by 

\An[i,M}\ 

M-*co M 

if the limits exists. If, moreover, the limit 

BD(A)= lim 

V ' M-L-^oo M - L 

exists (in which case it is certainly the same as D(A)), then it is called the Banach density 
of A(cf. gp.72]). 

The following lemma is routine. 

Lemma 1. Let (0 n )$£Li be a sequence of real numbers. Suppose N = [fj =1 Aj, where the 
union is disjoint. Let (#4 )n^i oe ^ e subsequence of (9 n ), consisting of those elements 9 n 
with n e Aj, 1 < j < r. 

1. If each (O^p) is asymptotically Fj -distributed for some distribution functions Fj, 1 < 
j < r, and D(A,) = dj, 1 < j < r, then (9 n ) is asymptotically F -distributed, where 
F = Ej=i djFj. 

2. If each {9^}) is asymptotically well Fj -distributed and BD(A,) = dj, then (9 n ) is 
asymptotically well F -distributed. 

Obviously, a general sequence on the line does not have to be asymptotically distributed 
according to some distribution function, but one would expect it of sufficiently "regular" 
bounded sequences. In our case, one might expect \n to be distributed according to some 
distribution function corresponding to a measure centered at about 1/C. However, this is 
not the case. In fact, the measure in question is supported on a finite interval, and is a 
convex combination of an absolutely continuous measure with decreasing density function 
on some interval [0, a] and the point mass 5 a at the right end a of that interval. 

We have defined E(N) indirectly by means of the following 

Problem 1. For each non- negative integer N, find non- negative integers 
n , m, . . . , riK-i, satisfying no + n% + . . . + nx-i = N, for which X^o* ^ ' n 1 IS minimal. 

It turns out that this problem is intimately related to a second optimization problem. 
Put Wi = 2£r, < i < K — 1, and let A denote the set of all positive odd multiples of the 
numbers 



A = {w , 3w , 5w , . . . , wi, 3ioi, 5t«i, . . . , w K - U 3w K ^, Bw K -i, ■ ■ ■} ■ 
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Here we treat A as a multi-set, or a sequence, in the sense that if some elements appear 
in this representation of A more than once (which occurs if some ratio Wi/wj is a rational 
number with odd numerator and denominator), then we consider A as having several copies 
of these numbers. 

Problem 2. For each non-negative integer N, minimize E™ = i $m, where 5i,5 2 , ■ ■ ■ ,5n 
range over all distinct iV-tuples in A. 

Note that, if an element appears several times in A, it is allowed to appear the same 
number of times in the sum as well. 

Let us demonstrate the equivalence of the two problems. Given the sum Y^Jq 1 Wi ■ nj, 
we may use the equality Wi ■ n? = Wi + 3wi + 5wi + . . . + (2rij — l)wi to see that any feasible 
value for the objective function of the first problem is a feasible value for the objective 
function of the second problem as well. On the other hand, solving Problem 2 is trivial. 
Namely, one minimizes the sum there simply by taking the N least elements of the set A. 
In particular, for each i, the multiples of Wi present in the optimal solution will be all odd 
multiples w iy 3wi, 5w iy ... up to some (2rij — l)u>j. Thus, the optimal solution of Problem 2 
yields the optimal solution of Problem 1 also. We note in passing that this discussion shows 
also that the minimum (for each of the problems) is obtained at a unique point unless A 
contains multiple elements. (However, we shall always refer to the optimal solution, even 
when there may be several.) 

A simple consequence of the above is 

Proposition 1. Let n = (n^fj^ 1 he the optimal solution of Problem 1 for some value of N. 
Then the optimal solution of Problem 1, with N + 1 instead of N, is n' = (n'^fs^ 1 , where 
n'j — n j + 1 for some < j < K — 1 and n\ = Hi for i ^ j. 

Remark. It is convenient to comment here on the effect a certain change in the original 
problem would make. One may consider the energies to be Wini(rii + 1) instead of Wiiif. 
This would change A to be the set of all even multiples of the w^s. Obviously, this would 
leave intact the equivalence of Problems 1 and 2. One can check that this would have also 
no effect on Theorems 1 and 2 infra. 

To formulate our main result we need a few definitions and notations. Real numbers 
6i,02, ■ ■ ■ ,0 r are independent over Q if, considered as vectors in the vector space R over 
the field Q, they are linearly independent. Equivalently, this is the case if the equality 
m\0\ + m 2 02 + . . . + m r r = for integer mi,m 2 , . . . ,m r implies m\ — m 2 — ■ ■ ■ — m r — 0. 
Considering the actual physical system (a collection of metallic grains), it is reasonable 
to assume that the capacitors Cj are random, so that generically they are independent 
over Q. Without loss of generality we may rearrange the K capacitors such that C = 
max <Kif_i Ci. It is also useful to divide all the capacitors by the largest one, so that the 
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scaled capacitors q = Cj/Co with 1 = Co > ci, C2 . . . , c^-i are dimensionless. Finally, set 
s = c + ci + . . . + c/r-i- 

Now we formulate our main results. 

Theorem 1. Suppose C ,Ci, . . . ,Ck-i are independent over Q. Tnen tie sequence 
(Xjv)jv=i j 's asymptotically F -distributed, where the distribution F is given by either of 
the following two representations: 



F(x) = {l-;E c iII( 1 



AT-1 ii"— 1 



i=0 



j=0 



2 Wjj 



x < 0, 

< x < 2w , 
2w < x, 



(4) 



= < 



0, 



x < 0, 



1 / r \ 1^1 

1-- £ (1^1 + 1)11^11(1-^)- !-^r > 0<x<2 Wo , (5) 

2w < 



SC{1,...,A--1} 



It is not immediately obvious from the formulas, but F has one discontinuity, namely at 
the point 2wq. The reason is that, as the elements of A are all odd multiples of the w^s, 
and as wq is the smallest of the w^s, it happens occasionally that there is no odd multiple 
of wi, . . . ,wk-i between two consecutive multiples of wq. The size of the atom at 2wq is 
- • n^=T 1 (l — This is easily explained intuitively. In fact, the "density" of odd multiples 
of Wi is Ci times the same density for multiples of wq. Hence the "probability" that an 
interval of the form [(2n — 1)wq, (2n+ 1)wq) does not contain an odd multiple of Wi is 1 — Cj. 
Assuming that the "events" of containing different tUj's are independent, we conclude that 
the proportion of multiples of Wq in A whose successors are also such is n^I 1 (l — q)- Since 
the proportion of multiples of Wq in A is 1, we arrive at the required expression for the size 
of the atom. 

Now we would like to study the asymptotic of the distances between consecutive elements 
of A as the number of capacitors grows. Obviously, as this happens, the distances become 
smaller. More precisely, on the average we have odd multiples of each Wj in each unit 
interval, and hence we have there y^^ 1 tt— = tt— elements of A altogether. Hence the 
average distance between consecutive elements is ^| fl . To understand the asymptotic of the 
gaps, it makes sense therefore to normalize them so as to have mean 1. Thus, we multiply 
the distances by and ask about the asymptotic behavior 

Theorem 2. Suppose the capacitances Cq,C±, . . . are chosen uniformly and independently 
in [0, 1]. For each K, let Fx denote the distribution corresponding to the normalized gaps 
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when taking into account the first K capacitors only. Then, with probability 1, the distri- 
butions Fk converge to an Exp(l) distribution function. 



Remark. As will be seen in the proof, we actually use much less to prove Theorem 2 than 
is required by the conditions of the theorem. Namely, we need the capacitances Cj to be 
linearly independent over Q, and that they do not form a fast diminishing sequence. 

It is worthwhile mentioning that this type of "Poissonian" asymptotic behavior of con- 
secutive gaps is typical. For example, this is the case for uniformly selected numbers in 
[0, 1], and is conjectured to be the case in other interesting well. (See, for example, 

|TT| and |L2| and the references there.) 

In the course of the proof, we shall make use of the notion of uniform distribution 
modulo 1 and a few basic results relating to it. (The reader is referred to Kuipers and 



Niederreiter ]10[ for more information.) A sequence (x n )^ =1 of real numbers is uniformly 
distributed modulo 1 if 

l{l<n<M:.< M < t)Ut _ 0<„<K1, 

M M-*oo ~ ~ 

where {t} is the fractional part of a real number t. In terms of Definition 1, (x n ) is uniformly 
distributed modulo 1 if and only if the sequence ({x„}) of fractional parts is F-distributed, 
where F is the distribution function of the uniform distribution on [0, 1]: 

{0, x < 0, 
x, < x < 1, 
1, X>1. 

The generalization of the notion of an asymptotically F-distributed sequence to that of 
an asymptotically well F-distributed sequence clearly carries over to our case. Instead of 
requiring only that the dispersion of large initial pieces of the sequence becomes more and 
more even, we require this to happen at arbitrary locations. This version is termed well- 
distribution. Thus, (i„)™ =1 is well- distributed modulo 1 if 

\{L<n<M:a<{x n }<b}\ 

— — > o — a, < a < o < 1 . 

M — L M-L^oo ~ ~ 

Both notions have multi-dimensional analogue. A sequence (x n )^L 1 in R s is uniformly 
distributed modulo 1 in R s if 

|{1 < n < N : a < {x„ } <b}\ 



H(bi-ai), 0<a<6<l 



where inequalities between vectors in R s are to be understood component-wise, 
(0, 0, . . . , 0) G R s , a = (ai, a 2 , . . . , a s ), and so forth. 
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Perhaps the most basic example of a sequence which is uniformly distributed modulo 1 
is (na)™ =1 , where a is an arbitrary irrational. In the multi-dimensional case, the sequence 
(nai, noLi-, ■ ■ ■ , na s ) is uniformly distributed modulo 1 in R s if and only if the numbers 
1, a%, «2, . . • , a s are linearly independent over Q. Moreover, in this case uniform distribution 
implies well-distribution (cf. |I0| , Example 1.6.1, Exercise 1.6.14]). 

Given a partition N = U*=i Aj and positive integers r 3 -, j = 1, we define the {rj) l - =l - 
inflation of the given partition as the partition of N obtained by inflating each element of 
each of the sets Aj into rj elements. More precisely, we construct sets Bj, j = 1, . . . ,1, as 
follows. For a positive integer i, let f(i) = j if % G Aj. Given any positive integer n, let m 
be defined by Yh^i f(i) < ni < Y%Li f(i) • Let n e Bj if m e Aj. The following lemma is 
routine. 

Lemma 2. In this setup: 

1. IfD(Aj) = dj, 1 < j < I, then D(Bj) = =P** 



2. If BD (A j) = dj, l<j<l, then BD (Bj) = . 

Proof of Theorem 1. Between any two consecutive odd multiples of Wq, there is at most 
one odd multiple of each Wj, 1 < j < K — 1. In fact, one easily verifies that, given a positive 
integer m, there is an odd multiple of Wj between (2m — l)u>o and (2m + l)u>o, namely there 
exists an integer n with 

(2m - l)w < (2n - l)wj < (2m + l)w , (6) 



if and only if 



mCj G 



l-Cj 1+Cj 



(modi). (7) 



Moreover, the relative position of (2n — l)wj within the interval [(2m — l)u>o, (2m + l)u>o) is 
the same, but in the opposite direction, as that of mcj (mod 1) within the interval (^ipS 
that is 

(2n - l)wj = a ■ (2m - l)w + (1 - a) ■ (2m + l)w , (0 < a < 1) , (8) 
if and only if 

mcj = (1 — a) ■ — — — + a ■ ° J (mod 1) . (9) 

Next we define a partition of N as follows. Write the elements of A in ascending order: 
A = {5i < 5 2 < 6 3 < . . .}. Given n e N, let S C {1, 2, . . . , K - 1} denote the set of all 
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those j's such that the unique interval of the form [(2m — 1)wq, (2m + 1)wq) containing S n 
contains an odd multiple of Wj. The set of all integers n giving rise in this way to any set 
S is denoted by B s . Consider the partition N = Usc{i,2,...,K-i} To prove the theorem 
using Lemma 1, we have to find the Banach densities of the sets Bs and the asymptotic 
distribution of the corresponding subsequences (Xn)neB s of Xn- 

The partition of N into sets of the form Bs is obtained as an inflation of a somewhat 
more straightforward partition. In fact, let S be any subset of {1,2, . . .,K — 1}. Denote 
by As the set of those positive integers n for which the interval [(2n — 1)wq, (2n + l)u>o) 
contains odd multiples of Wj for j G S and does not contain such multiples of the other w/s. 
Then N = \Jsc{i,...,K-i} As is a partition, and its (\S\ + l)sc{i,2,...,if-i}-inflation yields the 
partition N = \Jsc{i,2,...,k-i} b s- 

In view of the equivalence of (§) and (^), As is the set of those n's for which ncj G 
(~Y 2 ', ^y 2 - for j G S and ncj ^ : 4f I f° r j ^ S. By the conditions of the theorem, 

the numbers l,c±, . . . ,ck-i are linearly independent over Q, and hence the sequence c = 
(nci,nc2, ■ ■ ■ ,ncK-i)^=i * s well-distributed modulo 1 in H K ~ 1 . This means that 

D(A S ) = BD(A S ) = HciHil-d) . (10) 

■ies i^s 

Denote the right hand side of (|10D by ps- In view of the above and Lemma 2, this implies 

D(B S ) = BD(Bs) = + ^ — . (11) 

The denominator on the right hand side can be given a simpler form. In fact, let X iy i = 
1, 2, . . . , K — 1, be independent random variables with Xi ~ 5(1, q), and X = Y^f^i 1 -Xj- 
Then: 

(\T\ + l)p T = E(X + l) = l + Cl + ... + c K . x = s . (12) 

TC{1,2,...,A'-1} 

Hence: 

BD(£?5<) = M + VPs . (is) 

s 

Let 5* be an arbitrary fixed subset of {1, 2, . . . , K — 1}, say S = {1,2,...,/}, where 
< I < K — 1. If n G A s , then there exist odd integers a in , a 2n , . . . , a in such that aj n Wj G 
{(2n - l)w , (2n + l)w ). Put: 

v n = (ainWi, a 2n w 2 , a in wi) - (2n - l)w ■ (1, 1, . . . , 1) G [0, 2w ) 1 , n G A s . 

By the equivalence of (|) and ([5|), the sequence (v n ) ne A s is well-distributed modulo 2w in H l . 
Now each v n gives rise to / + 1 terms of (Xn)ne-B s > as follows. Let v$ < v^f> < . . . < v$ be 
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all coordinates of v n in ascending order. Set: 

The sequence (Xn)ns.z?s consists of of all coordinates of all vectors u n . Now we use the 
fact that if X\, X%, . . . , X r are independent random variables, distributed U(0,/i), and 
xw,xw,...,x^ are the corresponding order statistics, then each of the random vari- 
ables XW, X& - X^\ . . . , - X^-V, h - iW has the distribution function defined by 
G(x) = 1 — {x/h) r for < x < h (which follows as a special case from || p. 42, ex.23]). 
Consequently, for each 1 < j < I + 1, the sequence given by the jth coordinate of all 
vectors u n , n G As, is asymptotically well Gi-distributed, where Gi(x) = 1 — (x/2wo) 1 for 
< x < 2wq. Hence the sequence (Xn)neB s is asymptotically well Gi-distributed. Combined 
with (|T3|) , it proves (||). 

We shall indicate only briefly the proof of (|j), which is quite simpler. This time, we split 
(Xn) m to a union of subsequences (x$), < i < K — 1, by putting Xn in the sequence 
X$ if 5 n is a multiple of tOj. Clearly, the proportion of terms of (% n ) belonging to (Xn ) * s 
q/s. Next, consider the minimal odd multiples of all w^s which are larger than S n . The 
minimum of these K numbers is 8 n+ \. For each j ^ i, the distance from 5 n to the minimal 
odd multiple of Wj following 5 n is "distributed" U(0, 2vjj). (For i = it is also possible that 
the next term will be again a multiple of Wq.) The linear independence of the Cj's over Q 
implies that these K — 1 distances are (statistically) independent, so that their minimum is 
distributed according to the function G 2 (x) = 1 — n^o 1 (l — 2ur) on the interval [0, 2w ). 
These considerations can be formalized to prove (f|). This completes the proof. 

Remark. It is possible to shorten the proof by proving directly the equality of the right 
hand sides of (f|) and (^). In fact, it is easy to integrate both forms with respect to x; 
the equality of the resulting expressions follows easily from the binomial theorem. We have 
chosen the long way, as it is more instructive. 

Proof of Theorem 2. The distribution Fk is obtained from that in Theorem 1 by stretch- 
ing by the constant factor Hence: 



F K (x) 



0, x < 0, 

1 ( Cj x\ 

i-:E^II U-^r ' °< x<s ' (14) 



s i=Q j=o v s 
. 1, S < X 



Note that some of the values appearing on the right hand side depend on K implicitly. 
Namely, since w is assumed in Theorem 1 to be the least Wt, each time a is selected 
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which is larger than all the heretofore selected C/s, we have to rearrange the C/s, thus 
changing w and the c/s. We have to show that 



F K {x) — ► 1 - e~ x , x>0. (15) 

K — >oo 



Indeed, fix x > 0. Since 



s = c + ci + . . . + c^_i = — > Co + Ci + . . . + C K -i (16) 

and the Cj's are independent and uniformly distributed in [0, 1], we have 

S ^>oo. (17) 
Hence, with probability 1, for sufficiently large If we have 



1 K-l K-l / 

F K (x)=i-- y £c i n i 

S 4=0 3=0 \ 

3*/» 



2^ 



Thus, to prove (|T5|) we need to show that 



l^" 1 K ~} f c<X 



s £5 3=0 V s / 
Now, on the one hand, using the inequality 

1 - t < e~\ t e R , 

we have 

K ~ X f C X\ -^E?-o lfl 

III 1 -— )<e 8 < e- x+x ' s , 1 = 0,1,...,^-! 



3=0 



and therefore 



(18) 



E^.n i-- ?^. ^>o. as) 



7 E ^ n (l - -f) < \ E = e-*°" ^ e* . (20) 

6 i=0 3=0 \ 6 / 6 i=0 

On the other hand, as t — > we have 

e -(t+t 2 ) = ! _ (t + t 2 } + + 0(t 3 } = 1 _ t _ t_ + 0(t 3 } ) 

so that for all t in some sufficiently small neighborhood of 

e'^ <l-t. 

Consequently: 
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2 



j=0 



JJ ( i _ ^ ) > e ^ ' " > e -*-**l* . (21) 



Obviously, with probability 1, s grows linearly with X, namely for all sufficiently large K 
we have s > aK for a suitably chosen a > 0. (In fact, any a < | will do.) By (|2T|): 



j=0 



From @ and @ it follows that 



l^ 1 c ,x 



E^n i- 1 ? ^. ( 23 ) 



3'A 



which completes the proof. 
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